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Abstract 

We consider particle-antiparticle bound states in the scalar Yukawa (Wick-Cutkosky) model. 
The variational method in the Hamiltonian formalism of quantum field theory is employed. A 
reformulation of the model is studied, in which covariant Green's functions are used to solve 
for the mediating field in terms of the particle fields. A simple Fock-state variational ansatz is 
used to derive a relativistic equation for the particle-antiparticle states. This equation contains 
one-quantum-exchange and virtual-annihilation interactions. It is shown that analytic solutions 
of this equation can be obtained for the simplified case where only the virtual annihilation 
interaction is retained. More generally, numerical and perturbative solutions of the equation 
are obtained for the massive and massless-exchange cases. We compare our results with various 
Bethe-Salpeter-based calculations. 
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1. Introduction 

The variational method has been used sparingly in treating few-partiele bound and quasi-bound 
states in quantum field theory. Yet, it is in principle appealing, particularly for strongly coupled 
systems, because of its largely analytic and non-perturbative nature. The early papers that kindled 
interest in the variational method in QFT are those of Schiff [1], Coleman and Weinberg [2,3], 
Jackiw [4] Barnes and Ghandour [5], Stevenson [6], Tarrach [7] and others [8 ]. A brief review of 
this approach to few-body bound states in the Hamiltonian formalism of QFT up to 1996 is given 
inref. [9]. 

The variational method is, of course, only as good as the trial states that arc being employed. In 
the functional formulation (e.g. [5,8]) one is restricted largely to wave-functionals of the Gaussian 
type, due to the difficulty of handling analytically non-Gaussian functional integrals [10]. Another 
approach is to expand the trial state in a Fock-space basis. Indeed, the early work of Tamm [11] 
and Dancoff [12] is along these lines, though it was not variationally formulated. 

The Hamiltonian formalism of QFT can be expressed in terms of the QFTheoretic eigenvalue 
equation 

m) = E\^), (1) 

where H is the QFTheoretic hamiltonian operator, and E is the energy of the system under study. 
Such an equation is generally impossible to solve, except for some special models, particularly in 

1-|-1 (one spatial coordinate plus time), such as the Thirring [13 ]and Schwinger [14] models. In 
the variational approach within the Hamiltonian formalism one seeks approximate solutions to Eq. 
(1) by using the variational principle 

6{^t\H - E\^t) = 0, (2) 

where is a suitably chosen trial state containing adjustable features (parameters, functions). 
Some examples of the application of this method to bound states in scalar (j)^ theory are given in 
refs. [6c] and [15]. 

One of the difficulties of the Fock-state expansion approach for realistic models in 3-1-1 (such as 
the Yukawa model, QED, etc.) is that it results in an infinite system of coupled, multi-dimensional 
integral equations to be solved - an impossible task. Truncation, or other approximation schemes, 
undermine the strict variational nature of the approximation. This, then, puts into question the 
validity of results at strong coupling, precisely in the domain that one wishes to address in a 
non-perturbative approach. It has been pointed out recently [16] that various models can be 
reformulated in such a way that this difficulty can be circumvented. 

In this paper we shall implement the approach given in reference [16] to particle-antiparticle 
bound states in the scalar Yukawa model, in which scalar particles interact via a mediating real 
(massive or massless) scalar field. The treatment of two-body bound states in this model by means 
of the Bethe-Salpeter equation in the ladder approximation is known as the Wick-Cutkosky model 
[17,18]. In addition to these original solutions of Wick and Cutkosky, the scalar Yukawa model has 
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been used as a prototype QFT in many studies. It has been investigated quite extensively in various 
formalisms, such as the light-cone formulation [19-22], various Bcthc-Salpeter-based approaches 
[23-27], and others [28, 29, 30]. The work of Nieuwenhuis and Tjon [27b], in particular, gives a 
comparison of a number of quasipotential approximations. This makes the model appealing as a 
relatively simple test case. (We make no effort here to give an exhaustive survey of the literature 
on this model. Many relevant papers are cited in the references which we quote. A review of the 
Wick-Cutkosky model to 1988 is given by Nakanishi [31], while many mathematical details are 
given in the work of Silagadze [32].) 

In the scalar yukawa model massive scalar particles interact via a mediating real scalar field. It 
is based on the Lagrangian density {h = c = 1) 

L = {x)dy^{x)-rr?^*{x)(t>{x) (3) 

+ i^d^xix) dux{x) - ^f^^X^ix) - g(t)*{x)4){x) x{x) 

The mediating "chion" field can be massive ^ 0) or masslcss (i.e. ;U = 0, as in the original 
Wick, Cutkosky work [17, 18]). The coupling constant g has dimensions of (mass)2~T in N-l-1 
dimensionless. A slightly simpler model, in which ^ is real, is often considered. In that case there 
are only particles and no antiparticles. These models are closely related, since the forces among 
particles (and/or antiparticles) are only attractive (i.e. gravity-like, rather than electromagnetic- 
like). Of course, in the case where is real there is no particle-antiparticle annihilation. 



2. Reformulation of the model 



The fields (j) and x of the model (3) satisfy the Euler-Lagrange equations 

d^d,x{x) + ti\{x) = -g^*{x)^{x), (4) 

d''d,^{x) + mV(x) = -#(x)x(x), (5) 
and the conjugate of (5). Equation (4) has the formal solution 

X{x) = Xo{x) + j dx' D{x - x') p{x'), (6) 

where dx = d^xdt in + 1 dimensions, p{x) = (.x)(/>(x), Xo{x) satisfies the homogeneous (or 
free field) equation (cq. (4) with g = while D{x — x') is a covariant Green function (or chion 
propagator, in QFTheoretic language), such that 

(d''d^ + i?)D{x-x') = 5^+^{x-x'). (7) 

Equation (7) does not specify D{x—x') uniquely since, for example, any solution of the homogeneous 
equation can be added to it without invalidating (7). Boundary conditions based on physical 
considerations are used to pin down the form of D. 
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Substitution of the formal solution (6) into eq. (5) yields the equation 



d''d„(j){x) + m^(j){x) = -g4>{x)xoix) - gcpix) J dx'D{x - x')p{x'). (8) 

Equation (8) is derivable from the action principle 5 j dx C = 0, corresponding to the Lagrangian 
density 

jC = d''4>*{x)d,4>{x)-m^<t>*{x)<t>{x)-g(t>*{x)(t>{x)xo{x) (9) 
+ ]^ j dx' p{x)D{x - x')p{x'), 

provided that D{x — x') = D{x' — x). (We suppress the free chion part of the Lagrangian.) 

The QFTs based on (3) and (9) are equivalent in that, in conventional covariant perturbation 
theory, they lead to the same invariant matrix elements in various order of perturbation theory. 
The difference is that, in the formulation based on (9), the interaction term that contains the 
propagator leads to Fcynman diagrams involving virtual chions, while the term that contains xo 
correspond to diagrams that cannot be generated using the term with D(x — x'), such as those with 
external (physical) chion lines. 

The Hamiltonian density corresponding to the Lagrangian (9) is given by 



Hix) = n^ix) +n^{x) + ni,{x) + ni,ix), (lo) 



where 



and 



H^x) = (/)*(a;)0(x) + V(j)*{x) ■ V(pix) + (f)* {x)(t>{x) (11) 

Wx(^) = \xl + ^(Vxo)' + (12) 

Hi,{x)=g4>*{x)4>{x)x^{x), (13) 

Hi^{x) = -y / dx'(t)*{x)(t){x)D{x - x')<^*(x')<^(x'), (14) 



D(x-x')=f — e~^'='('^-'^') ^ fl5) 



where dk = d^~^^k and k ■ k = = k^k^. 
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To specify our notation, we quote the usual decomposition of the fields in A/" + 1 dimension: 

= J d^'q [(27r)^2c^,]-5 [^(q)e-*''-^ + St(q)e^'/-] (16) 

X{x) = j d^p[(27r)^2J)fc]-5[d(p)e-^P-^ + dt(p)e^P-^] (17) 

where Up = (p^ + m^) 2 , = (q'^ + /i^) 2 , q-x = q^Xy and q'^ = {q^ = Uq,q). The momentum-space 
operators A\ A, , B obey the usual commutation relations. The nonvanishing ones are 

[A(p),^t(q)] = [S(p),Bt(q)] =<5^(p-q) (18) 
[d(p),dt(q)]=<5^(p-q) (19) 

The Hamiltonian operator, H{t) = J d^x'H{x), of the QFTheory is expressed in terms of the 

particle and antiparticle creation and annihilation operators , A, B^ , B in the usual way. These 
operators are then commuted so that they stand in normal order in the Hamiltonian (we are not 
interested in vacuum-energy questions in this work). 

3. Ansatz for the particle-antiparticle system, variational equation and effective po- 
tentials. 

The simplest ansatz that can be chosen for a particle-antiparticle {(pfp) state is 

\H=J ciVd%F(pi,p2)At(pi)st(p2) I 0), (20) 

where F is an adjustable function. We use this trial state to evaluate the matrix elements needed 
to implement the variational principle (2), namely 



(V2|:i^</,--E:|V2) = J d''pid''p2F*{p^,p2)F{pup2)[ujp,+cOp,-E\, (21) 



and 



I d^pi d% dVl dV2 F*(p1, P^)F(P1, P2) 



9' 

8(27r) 

srJV / , „ „/ „/ \ +<^P2 -'^v'. ~'^vL )* ^ 



x6^{p, + P2- Pi - P^) e-^-^i-^-^^-^.-p^^- - (22) 

^iVp^Up^Upr^Up'^ 
1 1 1 



- (Pi + P2)^ - (p'l + ^2)^ - (p'l + Pi)^ - (P2 + ^'2)^ 
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We have normal-ordered the entire Hamiltonian since, at the present level of approximation (cf. 
the trial state (20)), this circumvents the need for mass renormalization which would otherwise 
arise in eq. (22). Also, in the Schrodinger picture we can take t = 0, and we do so henceforth. 

If we now specialize to the rest frame, where F(pi, P2) = /(pi)'^''^(pi +P2); then the variational 
principle (2) leads to the following momentum-space wave equation for the relative motion of the 
particle-antiparticle system: 



2ujp - E 



/(P) = 



4(27r: 



dV/(p') 



/«2 + (p' - p)^ - {^p 



(23) 



Note that the kernel (momentum-space potential) in this equation contains two terms. The first cor- 
responds to one-chion exchange and the second corresponds to virtual annihilation (this is perhaps 
more obvious from the four manifestly covariant terms in eq. (22)). 



In the nonrelativistic limit, 



EL 



<C 1, this equation reduces to 



PI 

m 



/(P) 



4(27r)^m2 



dV/(p') 



1 



/i2 + (p' - p) 



Am? 



(24) 



where e = E—2m. In coordinate space, equation (24) is just the usual time-independent Schrodinger 
equation for the relative motion of the particle-antiparticle system: 

- — V^V(r) + V{r)ij{r) = e^(r). (25) 

The potential V{r) is a sum of an attractive Yukawa potential (due to one-chion exchange) and an 
repulsive (if ^ < 2m) contact potential (due to virtual annihilation). In 3-1-1 dimensions these are, 
explicitly, 

V[r) = -a— + —^S^r), (26) 



where a 



4m2 — 

is the effective dimensionless coupling constant. 



IGvrm^ 

It is clear from Eq. (23) that the relativistic, momentum-space one-quantum exchange po- 
tential is always attractive, while the virtual annihilation potential is repulsive if the mass /j, of 
the mediating-field quantum is not too large, namely if /U < 2m. However, if /x > 2m then the 
annihilation potential becomes attractive at low momenta. 



4. Analytic solution of the variational particle-antiparticle equation with virtual an- 
nihilation interaction only. 

The variational particle-antiparticle equation (23) cannot be solved analytically. Of course, 
numerical solutions can be obtained, and these will be discussed in section 6. In addition, for 
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the massless-exchange case, analytic perturbative solutions can be worked out, and these will be 
presented in section 5. However, if only the annihilation interaction is kept (i.e. the chion-exchange 
interaction is turned off), exact analytic solutions of such a simplified particle-antiparticle equation 
can be obtained (for both bound and scattering states) . This is of interest as a solvable relativistic 
two-body equation, if for no other reason. Thus, if we neglect the first interaction term (the one- 
quantum-exchange term) in eq. (23) we obtain the rest-frame "annihilation-interaction" equation 



(27) 



where /o(p) is a "plane wave" solution of eq. (27) ( with g = 0) representing the particle and an- 
tiparticle incident on each other with energy E = 'iup^y Of course, /o(p) = for bound states. Be- 
cause the annihilation interaction is entirely repulsive if /x < 2m, we have only particle-antiparticle 
scattering solutions in that case. 

We shall discuss the solution of eq. (27) in some detail only in 3-1-1 dimensions. It is evident 
that the integral on the right-hand-side of eq. (27) vanishes, except in 5'-states, hence is sufHcient 
to write down the 5- wave component of eq. (27), namely 



fip) 



1 



S{p-po) 



where 



Prom equation (28) it follows that the S-wave phase shift r) is given by 



tan 77 = —Tram A- 



PO 



W2)2- 



Substitution of eq. (28) into eq. (29) yields the result 

1 



A 



{Att + aI{m,iJ,,po))' 



where 



roo 

I{m, ij,,pq) = m^V dp- 
Jo i 



p^ 



This principal value integral is, explicitly. 



(28) 
(29) 

(30) 

(31) 
(32) 



1 



2(g2 + 62) 1^ 
F=i= 2tanh-^ 



TT — 2 tan 



-1 



h 



+ 



IT q 



TT b 



vT + ^y b+l q{b+l) 



(33) 
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if < 6 < 1, where = 1 — (/x/2m)^, q = po/m, and 



+ 



+ l ' (g2 + 1)3/2 



(34) 



for the masslcss exchange case, /x = (i.e. 6=1). With A as given in eq. (34), the tangent of the 
S-wave phase shift becomes 



tan rj = — 



ira m po 



(ul^-{li/2Y) (47r + «/(m,/x,po))' 



47r 



(35) 



from which the elastic particle-antiparticle scattering cross section a = — sin r] is readily calcu- 
lated. 

For the masslcss chion exchange case (^u = 0), the cross section, in units of vr/m^, starts from 



a value of 



G-ln- 



(16 + [l + 4/7r]Q;)2 
a2 



aX q = PQ/m = then decreases monotonically with increasing q to 



the asymptotic form — as g — oo. Note that the maximum value of the cross section (which 

occurs at = for all a in this massless-exchange case) increases uniformly from zero at a = to 
an asymptotic value of 647r^/(7r + 4)^ = 12.3848 as a ^ oo. 

When the mediating field quanta are massive then, for given a, the cross-section, as a function 
of collision energy, behaves qualitatively in a similar way to the massless case if /x/m < 1. However 
the shape of the cross section changes as /x increases towards 2m, in that it initially increases with 
the collision energy, reaches a maximum and then decreases towards zero as 
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1 1 a' 

4^ ^ 16 (6 + l)q 



3 1 

7+0(-8)- 



(36) 



We shall not delve into a detailed discussion of the behaviour of a{q) for various a and h since the 
specific results can always be evaluated using the given analytic formulae (33 - 35). 

It is of interest to note that in the non-relativistic limit, i.e. p <C m, eq. (27) becomes. 



e 

m 



/(P) = /o(p) + 



4(27r)^ m2(^2 



(37) 



where e = E — 2m. The resulting / integral (cf. eq. (32)) diverges in N = 2 and 3 dimensions, 
thus resulting in a vanishing phase shift and cross section. This just reflects the "trivial" nature 
of the scattering by a rcpiilsive delta function potential in non-relativistic (Schrodinger) theory in 
A'^ > 1 dimensions ("trivial" in the sense that the ^-matrix is unity). 
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For very massive mediating fields, fi > 2m, the annihilation interaction becomes attractive at 
low momenta, and this leads to binding of the particle-antiparticic system if the coupling constant 
is large enough. The energy eigenvalue condition is (from eq. (27) with /o = 0) 

''d^'pj:^ r^. (38) 



4(27r)^J ^(2a;p-£;)u;2(^2_4^2) 
In iV = 3 dimensions equation (38) yields the result 



27r 2 \ b . ^ f b \ , / ^4 

- tanh"^ , + tan"^ ' 



e 



2 



a 4 + 4b^-e^[^/lT¥ \VlT¥ 

b"^ e 1 TT 1 1 / 4 . 

+ , .2 + 7;tT-^- 1 - -7^=^ , (39) 



1 + 62 ^/izr^ j ' 2 1 + 62 e V V4 



where, now, 6^ = (///2m)^ — 1 and £ = £'/m. A plot of E{a) for a representative value of 6 = 1, 
i.e. fi = 2\/2m, is given in Figure 1 (actually, a{E) is plotted). Note that binding does not set in 
until the coupling constant a exceeds a minimum value ag. Thereafter the energy (i.e. particle- 
antiparticle mass) decreases monotonically with increasing a to an asymptotic value E/m = Smin, 
where both ao and emin vary with /i/m. The general expression for ao is 



ao = 27r 



1 1 , 1 / 6 \ TT 1 

tanh"^ -== + -- 



2 6V62TT KV^TTJ 4 62 + 1 



TT + 2 3(7r + 2)2 



The minimum value of a at which binding sets in (vis. ao = = 4.88812) occurs when 6 = 
(i.e. when /x just passes 2m), but the binding is very weak, since emin is barely below 2 in that 
case. As 6 increases, so does ao and so does the binding energy. For a near ao the behaviour of 
E{a) is of the form 




+ Oi{a-aof). (41) 

Once again, it is of interest to note that the annihilation interaction, in A'' = 2,3 spatial 
dimensions, can support bound states only in the relativistic formulation. This is because, had we 
started from the non-relativistic form of eq. (37) (with /q = for bound states) , then the eigenvalue 
equation corresponding to equation (38) would be 



1- »^ 



['fp-j^ (42) 

J (p'^ m — e) 



4(27r)^ m2(/i2 _ 4^2) J ^ (^2/^ _ ' 

The integral (42) converges only for A" = 1, whereas in the relativistic case the integral (38) 
converges for N = 1,2,3. In other words the short-range annihilation interaction does not sup- 
port bound states non-relativistically in 3 or even 2 spatial dimensions (the interaction is an at- 
tractive delta-function potential in the non-relativistic limit), but it does support bound states 
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relativistically. This implies that the virtual annihilation interaction strengthens (relative to the 
non-relativistic delta function potential) if relativity is taken into account. 



5. Perturbative results in 3+1 dimensions for the particle-antiparticle binding energy 
in the massless-exchange case. 

The relativistic two-particle equation (23) can be reduced to radial form by setting 

/(P) = f{p)YUp), (43) 

where p = |p| and y^m(p) are the usual spherical harmonics, and carrying out the angular integra- 
tion. The result, in iV = 3 dimensions, is 



with 



where 



2lo„ - E 



z = 



f(p)=- r dp'?-f{p')K,{p',p), 

Jo P 



m 



Qe{z) 



vr 



^-^PP'^eo 



2pp' 



(44) 



(45) 



(46) 



a = and Qe(z) is the Legendrc function of the second kind. This equation is similar to 

that derived by Di Leo and Darewych using a variational-perturbative approach [33] . That previous 
result corresponds to eq. (45) without the virtual-annihilation interaction (the second term in eq. 
(45)), and with Up = ujp' in the first, one-chion exchange term of eq. (45). 

Since the solutions of equation (44) in the non-relativistic limit (without the virtual annihilation 
interaction) are the well known hydrogenic wave- functions in momentum-space [34], we can use 
them to obtain perturbative expressions to the particle-antiparticle mass (rest-energy). The result 
is 



1 



1 



2m — -ma — ^ ma 



+-ma 



16 
4 



3 
4^ 



{2i + l)n3 n 



{2e + l)n3 

^ ^ + AEanni + 0(a')- 



(47) 



The terms on the right are the rest energy, the non-relativistic Balmer term, the O(a^) correc- 
tion to the kinetic energy, one-chion exchange interaction and the virtual annihilation interaction, 
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respectively. The correction due to the annihilation interaction is 



^Eanni = -mM / dp dp' f {p)f ip')-^p'' p'^ Seo 

= lma^^Seo + 0{a^) (48) 

The result (47) agrees with the earlier work [33], except that the annihilation correction (48) is 
new. The virtual annihilation correction is reminiscent of what is obtained for triplet i = states 

of positronium, where one obtains —ma^^6ioSsi, with a = e^/hc in that case. Note that the 

4 n-^ 

"retardation term", (uop — oopi)"^ of Eqs. (23) and (46), has no effect at 0{< a^). 

As is well known, the massless Wick-Cutkosky model has been solved in the ladder approxima- 
tion of the Bethe-Salpeter formulation [17, 18], as well as the light-cone ladder approximation [19, 
20]. The expansion of these solutions in powers of a is found to be [19] 

E/m = 2-^- + 0(a3) (49) 

4n^ im^ 

This is quite different from our result (47) beyond the 0{a'^) Balmer term. The unusual 0{a^ In a) 
terms are an artefact of the ladder Bcthe-Salpctcr formulation, and allegedly do not arise if crosscd- 
ladder diagrams are included [35]. A discussion of the origin of the a'^lna term is given by A. 
Amghar and B. Desplanques [36]. In any case, our results contain no such terms, and are much 
more like the corresponding results for positronium in this respect (i.e. that the lowest order 
relativistic corrections to the Balmer result are 0{a^)). 



6. Numerical solution for the ju/m = and /x/m = 0.15 cases. 



We have solved equation (44) approximately for the ground state in = 3 spatial dimensions, 
ng the variational method with the ti 

parameter, determined by minimizing 



using the variational method with the trial wave function f{p) = ^ ^ .^iji^n ' ^ adjustable 



E 



poo /y rco POO 

/ dp2upp''f\p)-- / dp'dpp'pf{p')fip)Ke{p',p) 
Jo vr Jo Jo 



rco 



'0 



dpp'fip), (50) 



with respect to b, for various n and given a. Although these variational results are only approxi- 
mations to the "exact" (i.e. numerical) solutions of eq. (44), they are in fact reasonably close to 
the numerical ones for the entire range of values of a considered, as shown in a previous study 
[33]. A list of ground-state values of E/m in the massive-exchange case for the present model is 
given in Table 1, for n = 2,2.5,3,3.5. Generally, the n = 2 values are lowest for a < 0.5, where 
relativistic effects are not so pronounced, but the n = 2.5 and n = 3 values are lower at strong 
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coupling. However, the various values look quite similar on a graph, and in the figures we shall plot 
curves corresponding to a single value of n only, as this will be sufficient for comparison purposes. 

An advantage of using the variational solution is that we have an analytic representation of the 
wave-function, and so can examine its behaviour as the coupling constant a changes. Thus, we note 
that the values of the "inverse Bohr radius" parameter b, for any given n, increase monotonically 
with increasing a. Also, for any given value of a, the parameter b increases as n increases. 

A plot oi E{a, fi = 0) for the ground state obtained in this way is shown in Figure 2. We plot two 
versions of our results, namely with and without the virtual annihilation interaction (second term 
of eq. (45)) included. Note that the affect of virtual annihilation is substantial, and it increases 
with increasing a. 

As mentioned previously, the scalar Yukawa (or Wick-Cutkosky) model has been studied by 
many authors in various formalisms and approximations. It is therefore of interest to compare 
some of them to our results. Thus Figure 2 also contains plots of the classic solutions of Wick 
[17] and Cutkosky [18] of the Bethe-Salpeter equation in the ladder approximation, as well as the 
analogous light-cone calculations of Ji and Furnstahl [20]. We also plot the results of Di Leo and 
Darewych [33], which correspond to the present results with ujp = LOpi (i.e. no retardation). None 
of these results contain the virtual annihilation interaction (which is repulsive) and so they should 
be compared with the present results without virtual annihilation. 

It is evident from Fig. 2 that our variational results predict stronger binding that either of 
the ladder Bethe-salpeter results, or the Di Leo results. (The difference between our no-virtual- 
annihilation results and the Di Leo results show the effect of the retardation term {ujp — uJp'Y in 
the potential of eq. (23).) Indeed, our results predict stronger binding than ladder B-S, even if 
we include the repulsive virtual annihilation interaction. Numerical values corresponding to Figure 
2 are listed in Table 2 for n = 2, along with the wave-function parameter b. Table 3 is a list of 
corresponding results for n = 3. 

It has been shown recently [37] that exact two-body eigenstates can be written down for the 
QFTheoretic Hamiltonian if^ + Hi^ (cf. eqs. (10, 14)), that is, for the present model without 
free chions, provided that an "empty" vacuum state |0), annihilated by both positive and nega- 
tive energy components of the field operator ^(x), is used. The use of such an empty vacuum 
state results in a relativistic two-body scalar equation that has both positive and negative energy 
solutions. Such negative-energy solutions do not arise (and should not arise) in the conventional 
QFT treatment that uses a Dirac "filled-negative-energy-sea" vacuum state, including the present 
work. Nevertheless the two-body equation obtained in [37] can be solved analytically [37, 38] in 
the massless-mediating-field case, and the positive energy, ~ 2m' - like solution is 



E = m 



i 




(51) 



This result is quite different from that of the present treatment, or the ladder Bethe-Salpeter 
calculations. For one thing there is a rather low critical value of the coupling constant, namely 
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CKc = n (oc = 1 for the ground state) beyond which the two-body energy (rest mass) ceases to be 
real {E/m = = 1.414... when a = a^)- This is similar to what happens for one-body Klein- 
Gordon or Dirac equations in a Coulomb potential. We also plot the result (51) in Fig. 2, and note 
that in the domain < a < 1 it predicts stronger binding than any of the other results shown in 
the figure. 

(An aside: The formula (51), and its mi ^ m2 generalization derived in ref. [38], was obtained 
previously by Todorov [25] using a quasipotential approach. Todorov's article [25], which came to 
our attention recently, also contains a useful historical overview of earlier work on the relativistic 
two-body problem in QFT, including references to various "rediscovered" formulations and results.) 

Figure 3 is a plot of the particle-antiparticlc ground state energy in = 3 spatial dimensions for 
the massive exchange case, with fi/m = 0.15. Once again, we plot the solutions of our variationally 
derived equation (44) with and without virtual annihilation included. In addition, we plot some 
Bethe-Salpeter based quasipotential results that are given in a study of the 4>'^x model by Nieuwen- 
huis and Tjon [27b], together with their numerical results obtained using a Feynman-Schwinger 
formulation. The various quasipotential results plotted in Figure 3 are explained in the paper of 
Nieuwenhuis and Tjon [27b], and this will not be repeated here. None of the calculations, save 
ours, contain the virtual annihilation interaction, so they should be compared to the version of our 
results without virtual annihilation. Some numerical values corresponding to the curves plotted in 
Figure 3 are listed in Table 4. 

It is evident from Figure 3 that the ladder Bethe-Salpeter (BS) calculation predicts the weakest 
binding, while the Feynman-Schwinger (FS) results of Nieuwenhuis and Tjon predict the strongest 
binding. (The FS calculations arc nonpcrturbativc, but contain no loop effects.) The various 
quasipotential results arc distributed between the ladder BS and the FS values, in the following 
order: Blankenbecler-Sugar-Logunov-Tavkhelidze [39], Gross (with retardation) [40] and "equal- 
time" [41]. All these quasipotential results are taken from Figure 1 of Nieuwenhuis and Tjon [27b]. 
Our present results (without virtual annihilation) fall above most of the quasipotential ones, and 
lie closest to the Blankenbeclcr-Sugar curve. We also include a plot of the two-body Klein-Gordon 
Feshbach-Villars formalism results of ref. [37], which lie very close to the Gross results. We can 
only speculate to what extent this is a coincidence (the Gross and FV formalism equations are not 
obviously similar). 

Thcussl and Desplanques [42] have recently calculated the two-body energy being studied here. 
They use the Bethe-Salpeter equation, but include crossed-ladder effects in an approximate per- 
turbative way. Their results fall very close to the Gross quasipotential values (and to those of the 
Feshbach-Villars formulation) in the domain a < 1. Therefore, we do not plot them in figure 3. The 
Theussl-Desplanques work underscores the inaccuracy of the ladder Bethe-Salpeter approximation 
and the importance of including crossed-ladder effects. 

The fact that the present variational results show considerably weaker binding than the FS 
results of Nieuwenhuis and Tjon is perhaps not surprising, since the present calculation contains 
no "crossed-ladder" effects (the simple variational ansatz (20) is incapable of incorporating such 
effects). However, our results show considerably stronger binding than the ladder BS, even though 
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the latter uses essentially the same kernel (one-chion exchange) as the present calculation. 

The quasipotential results differ from our variationally derived values in several respects. For 
one thing, the quasipotential equations are all different and somewhat ad-hoc (though physically 
motivated) one-time modifications of the BS equation, whereas our equations are obtained in a 
completely "ab-initio" way, and are limited primarily by the simple choice of ansatz (20) that we 
have made in this work. The quasipotential equations, in the unequal mass case, have the Klein- 
Gordon (KG) equation as their one-body limit. In this sense they arc more like the two-body KG 
Feshbach-Villars results of refs. [37, 38], which also have the KG eq. as their one-body limit. The 
present variational approach leads to an equation, which, in the unequal mass case, does not have 
the KG one-body limit. 

The stipulation that the "correct" one-body limit of a relativistic two-scalar-particle equation 
should be the KG equation is often made in the literature. However, this is a curious measure 
of correctness, since the KG equation has negative-energy solutions, which should not arise in 
a conventional QFT treatment that uses a Dirac "filled-negative-energy sea" vacuum. Indeed, 
negative energy solutions cannot, and should not, arise in the present calculations (we use the 
conventional Dirac vacuum), and so the unequal- mass counterpart of eq. (27) does not (and should 
not) have the KG equation as its one-body limit. On the other hand, the two-body KG Feshbach- 
Villars formalism equation (given in ref. [38] for the unequal mass case) is obtained using an 
"empty" vacuum, and so it does have negative-energy solutions, and also the KG one-body limit. 
Its binding energy predictions are similar to those of the quasipotential equations, even though 
it contains no retardation, or manifest crossed-ladder effects. The one-body limit equation in the 
present formalism has only positive-energy solutions, like the two-body equation (27). 

Of course, we do not claim that our present results arc "better" than the quasipotential equation 
ones (in the sense that they are closer to the unknown exact results for this model). The trial state 
(20) that we are using here is too simple to make any such claim. However, our approach is strictly 
variational, with nothing "put in by hand", and can be systematically improved by improving on 
the trial state (20) (which we are in the process of doing). It may be that the numerical FS results 
of Nieuwenhuis and Tjon are the most accurate binding-energy results available to date for the 
scalar Yukawa model, so it will be interesting to see how improved variational results will compare 
with this benchmark. 

We should point out that the trial state (20) that has been used in this work is insensitive 
to the if/j term of the Hamiltonian (cf. Eq. (13)). Thus, it is suitable for describing stable 
particle-antiparticle states only, without explicit annihilation of the particle-antiparticle system, or 
decay of the excited bound states, with the emission of physical chions. It is possible to include 
such processes in the present formalism in various ways, such as perturbatively, or by suitable 
modification of the ansatz (20) (see, for example, refs. [33], [43]). However, we do not consider 
such processes in this paper. 
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7. Concluding remarks. 



We have applied the variational method to the study of particlc-antiparticle bound states in 
the scalar Yukawa (Wick-Cutkosky) model (scalar particles interacting via a massive or massless 
mediating scalar "chion" field). We have used a reformulated version of this theory in which a 
covariant Green function is used to eliminate the chion field partially, so that the chion propagator 
appears directly in the QFTheoretic hamiltonian. 

A simple Fock-space trial function is used in the variational method. It leads to a relativistic 
particle-antiparticle momentum-space equation with the covariant one-chion exchange and virtual 
annihilation Fcynman amplitudes appearing in the kernel (momentum-space potential) of the equa- 
tion. The virtual annihilation interaction is repulsive, except if the mcdiating-ficld quantum (the 
chion) is very heavy (more massive than the combined rest mass of the particle and antiparticle) 
whereupon it becomes attractive at sufficiently large momenta. This particle-antiparticle equation 
has no negative energy solutions, i.e. it is free of any negative-energy "pathologies" . 

When the one-chion exchange interaction is turned off, we find that the resulting model theory 
(with a purely virtual-annihilation interaction) is analytically solvable for scattering states, and also 
for bound states when the chions are very massive. The virtual annihilation interaction reduces 
to a delta-function (contact) potential in the nonrelativistic limit. It supports no bound states 
in A'^ = 3 spatial dimensions and the S'-matrix is unity in this limit. However, if the relativistic 
equations are used, we find that bound states are possible and the scattering is not trivial. This 
analytically solvable relativistic model is instructive in understanding the effects of a relativistic 
generalization of the delta-function potential, such as occurs in the virtual annihilation interaction. 

In the general case, with both one-chion exchange and virtual annihilation interactions included, 
the relativistic particle-antiparticle equation cannot be solved analytically (at least, we do not know 
how to do so). However, analytic expressions for the energy (rest mass) of the bound particle- 
antiparticle system can be obtained as an expansion in the effective dimensionless coupling constant 
a for the case of a massless mediating field. We give such an expression for the ground and arbitrary 
excited states of the system to 0(q!^) inclusive. We find that the lowest-order relativistic corrections 
to the Balmer formula are 0(ol^)^ much like for positronium, and quite unlike the predictions of 
the Bethe-Salpeter equation in the ladder approximation, which include unusual o? In a terms. 

Lastly, we calculated the particle-antiparticle bound-state energy for arbitrary a in the ground 
state. This was done using a variational approximation rather than numerical integration of the two- 
body equation, as the results are not much different and the variational method allows one to exhibit 
the behaviour of the wave-function more transparently for various values of a. We study the case 
where the exchanged quantum has mass /j, = 0, and find that our results predict stronger binding 
than ladder Bethe-Salpcter approximations (Wick-Cutkosky or light-front solutions). For the case 
of massive chion exchange, with /x = 0.15, we find analogous behaviour. However, our results show 
weaker binding than most of the quasipotential reductions of the Bethe-Salpeter equations. 

The present approach has several attractive features. Firstly, it leads to equations with no 
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negative-energy or mixed-energy solutions, such as arise in many other formulations (though such 
negative-energy and mixed-energy solutions usually are not discussed). Secondly, the results are 
strictly variational, with no pcrturbative approximations. Thus the results are applicable at strong 
coupling, at least in principle (we hasten to add that variational results are only as good as the 
trial states employed in their use). Unlike the quasipotential reductions of the Bethe-Salpeter 
equation, our results are rigorous in the sense that nothing is put in by hand. Thirdly, the method 
is amenable to systematic improvement by improving the variational trial state. Lastly, the method 
is straightforwardly generalizable to relativistic three or more particle systems, though, of course, 
one is then faced with the usual complexity of a relativistic many-body problem [44] . 

We thank B. Dcsplanques for sending us preprints of recent work, and for useful conversation. 
The support of the Natural Sciences and Engineering Research Council of Canada for this work is 
gratefully acknowledged. 
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Figure captions 



Figure 1 

Two-particle bound-state energy E (eq. 39) due to virtual-annihilation interaction only, for 
various values of the coupling constant a, at a mediating-quantum mass of /j, = 1\f2m (i.e. 6=1). 
Binding sets in for a > ao = 8.9203, and E/m ^ emin = 1-5058 for nlm = 2\/2. 

Figure 2 

Particle-antiparticle bound state mass Elm for various values of the coupling constant a, mass- 
less quantum exchange. Curves, from highest to lowest: ligh-front ladder Bethc-Salpeter [Ji and 
Funstahl, 20a]; ladder Bethe-Salpeter [Wick and Cutkosky, 17, 18]; equations (44) without virtual 
annihilation or retardation [Di Leo and Darewych, 33]; equation (44) with virtual annihilation 
and retardation (n = 2); equation (44) with retardation but without virtual annihilation (n = 2); 
Feshbach-Villars formalism [Darewych, 37]. Note that all these results, save one, do not contain 
the virtual-annihilation interaction. 

Figure 3 

Particle-antiparticle bound state mass Ejm for various values of the coupling constant a, mas- 
sive quantum exchange: ^u/m = 0.15. Curves, from highest to lowest: ladder Bethe-Salpeter [Wick 
and Cutkosky, 17, 18]; present results of equation (44) with virtual annihilation {n = 2); equa- 
tion (44) without virtual annihilation (n = 2); crosses: Blankenbecler-Sugar [27b, 39]; Feshbach- 
Villars formalism [Darewych, 37]; diamonds: Gross [27b, 40]; boxes: equal-time [27b, 41]; circles: 
Feynam-Schwingcr formalism [27b]. Note that all these results, save one, do not contain the virtual- 
annihilation interaction. 
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Table 1. 



Ground state particle-antiparticle mass (rest energy) E^jra obtained from eq. (50) using /(p) = 
LOpI {p^ + 6^)" for various values of the coupling constant a = ^^^^i and variational parameter n. 
Massive exchange case, ji/m = 0.15. Virtual annihilation is included here. The number in brackets 
below each energy is the corresponding value of 6/m. 
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Table 2. 



Value oi E2/m (ground state) for ii/m = in 3 + 1, eq. (50) with n = 2. 
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Table 3. 



Value oi E2/m (ground state) for ii/m = in 3 + 1, eq. (50) with n = 3. 
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Table 4. 



Value of E2/m (ground state) for /i/m = 0.15 in 3 + 1. Results given in columns 2,3,4 and 6 
were read off fig. 1 of [27b], hence the accuracy of the last figure is questionable. 



a Niewenhuis Equal Gross eq. Darewych Blankenbecler Di Leo Present results 
and Tjon Time (with retard.) [37] and Sugar [33] (n = 2) no annih. 



0.3 








1.999536 




2.0 


2.0 


0.4 


1.99 






1.99534 


1.996 


1.996946 


1.996582 


0.5 


1.98 






1.98630 


1.989 


1.990833 


1.989742 


0.6 


1.96 


1.966 


1.969 


1.97176 


1.979 


1.982148 


1.979868 


0.7 


1.91 


1.941 


1.948 


1.95081 


1.965 


1.971207 


1.967249 


0.8 


1.85 


1.907 


1.919 


1.92199 


1.952 


1.985276 


1.952143 


0.9 


1.77 


1.861 


1.880 


1.88282 


1.933 


1.943577 


1.934776 


1.0 








1.82847 




1.927296 


1.915346 


2.0 












1.703450 


1.637352 


3.0 












1.412621 


1.262094 


4.0 












1.084915 


0.830372 


5.0 












0.733877 


0.361942 



23 



10000 r 1 i i 1 i i i 

: b(x) 



1000 r 



100 r 



10 r 




■\ I I I I I I I I I 

1.4 1.5 1.6 1.7 1.8 1.9 2 2.1 2.2 

E/m 



